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THE  INTERACTION  AND  DAMAGE  EFFECTIVENESS  OF  TWO 
EXPANDING  NON-SPHERICAL  EXPLOSION  BUBBLES 


1 .  INTRODUCTION 


The  underwater  explosion  of  two  simultaneously  or  sequentially 

initiated  separated  charges  produces  complex  shock  wave  and  bubble 

phenomena.  The  interaction  of  the  two  shock  waves  has  been  studied 
1  2 

in  the  past  ’  .  This  report  is  a  continuation  of  a  study  of  the 
bubble  phenomena  of  separated  charges. 


With  two  separated  charges,  two  distinct  bubbles  are  formed 


following  the  detonation.  The  bubbles  will  deform  from  a  spherical 
shape  as  a  result  of  the  hydrodynamic  interaction.  The  authors 
have  shown  that  the  centers  of  two  spherical  explosion  bubbles 
remain  essentially  stationary  during  the  expansion  phase  but  later 
attract  strongly  as  the  bubble  contracts.  It  is  assumed  that 
explosion  bubbles  which  are  not  constrained  to  a  spherical  shape 
will  behave  in  the  same  manner.  The  analysis  deals  with  explosion 
bubbles  that  pulsate  in  phase.  Only  the  expansion  phase  will  be 
considered.  The  effect  of  buoyancy  upon  the  bubbles  is  neglected. 

Our  study  utilizes  the  methods  of  potential  flow  theory  to 
describe  the  expansion  of  the  bubbles  and  their  shapes  as  a 


N.  L.  Coleburn  and  L.  A,  Roslund,  "Collision  of  Spherical  Shock 
Waves  in  Water,"  NOLTR  68-110,  (1968) 

2 

R.  J.  beeger  and  H,  Polachek,  "On  Shock  Wave  Phenomena:  Waterlike 
Substances,"  Journal  of  Applied  Physics,  22,  640  (1951) 

R.  G.  Leamon,  C.  G.  Simpson,  and  H.  G.  Snay,  "The  Interaction  of 
Two  Spherical  Explosion  Bubbles,"  NOLTR  73-101,  (1973) 
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function  of  time.  Lagranglan  coordinates  are  used  to  follow 
individual  particles  on  the  bubble-water  Interface  and  for  the 
determination  of  the  shape  of  the  bubble.  This  problem  has  been 
studied  elsewhere  using  the  far  more  sophisticated  MACYL  computer 
technique.  The  results  of  our  approximation  compare  favorably  with 
those  obtained  with  the  MACYL  analysis  that  uses  finite  difference 
techniques.  The  potential  flow  calculations  are  simple  and  require 
only  a  small  portion  of  the  machine  time  needed  for  the  finite 
difference  method. 


2 .  METHOD 

In  the  method  described  below,  a  solution  of  the  appropriate 
hydrodynamic  equation — Laplace's  equation — is  used.  Such  a  solution 
specifies  the  particle  velocities  in  the  water.  The  initial  shape 
of  each  bubble  and  the  equation  of  state  for  the  water  and  that  for 
the  gaseous  contents  of  the  bubbles  are  those  of  the  classic  bubble 
theory,  namely  that  of  an  incompressible  medium  for  the  water 
surrounding  the  bubbles  and  that  of  an  ideal  gas  for  the  contents 
of  the  bubbles. 

The  bubble  Interface*  is  assumed  to  consist  of  the  same 
particles  of  water  throughout  the  pulsation.  If  the  particle 
velocities  in  the  water  are  known,  and  if  the  initial  Interface  is 
specified,  then  the  shape  of  the  bubbles  can  be  calculated  as  a 
function  of  time  by  following  the  motion  of  the  Interface  particles. 

This  approach  allows  the  bubbles  to  combine  or  coalesce. 

If - 

For  the  purpose  of  this  report,  the  term  Interface  shall  refer 
to  the  bubble-water  Interface  regardless  of  whether  or  not  the  two 
bubbles  are  coalescent. 
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The  Velocity  Potential.  For  the  purpose  of  describing  the 
motion  of  the  bubbles,  the  flow  field  In  the  water  Is  taken  to  be 
Irrotational  and  Incompressible.  These  assumptions  form  the  basis 
of  the  well-established  bubble  theory,^’^  the  classic  bubble  theory. 
Under  these  conditions,  a  velocity  potential,  (|),  exists  which 
describes  the  flow  field  and  which  satisfies  Laplace’s  equation 

^^<}>  =  0  .  (2.1) 


Any  velocity  potential  which  satisfies  this  equation  Is  a 
possible  solution  to  our  problem.  One  class  of  functions  which 
satisfy  Equation  (2.1)  are  called  singular  solutions,  so  called 
because  (p  becomes  infinite  at  the  origin  of  these  solutions.  These 
include  sources,  dipoles,  quadrupoles,  etc.  The  simplest  of  these 
solutions,  the  source.  Is  written  as 


where  Q(t)  Is  an  arbitrary  function  of  time  and  r  Is  the  distance 
from  the  origin  to  any  point  In  the  flow  field.  In  fact  the  source 
Is  the  complete  solution  of  a  single  pulsating  non-migrating 
spherical  bubble.  In  this  case  It  takes  the  form 


B.  Friedman,  "Theory  of  Underwater  Explosion  Bubbles,"  Institute 
of  Mathematics  and  Mechanics,  New  York  Universltj'',  Report  No. 
IMM-NYU_l66,  (1947).  Reprinted  In  Underwater  Explosion  Research. 
Volume  II,  The  Gas  Globe,  Office  of  Naval  Resear'^,“p;  329-37B,  (1950) 

Snay  and  E.  A.  Christian,  "Underwater  Explosion  Phenomena? 

The  Parameters  of  Non-Migrating  Bubbles,"  NAVORD  2937,  (1952) 
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(2.3) 


where  A  Is  the  radius  of  the  spherical  bubble  and  A,  its  time 
derivative  which  is  the  radial  velocity  of  the  bubble  interface. 

We  use  two  equal  source  potentials  fixed  in  the  flow  field  to 
represent  the  motion  due  to  the  two  bubbles.  These  bubbles  are 
taken  as  being  formed  by  equal,  simultaneous  underwater  explosions. 
The  effect  of  buoyancy  on  the  bubbles  is  neglected.  This  effect  is 
known  to  be  small  for  a  single  bubble  during  the  expansion  phase^. 
In  addition  there  are  compelling  reasons  for  omitting  the  effect  of 
buoyancy*.  Using  the  axlsymmetric  coordinate  system  shown  in 
Figure  1,  the  form  of  the  velocity  potential  becomes 


Where  B  is  the  total  separation  distance  between  the  two  charges. 
Physically,  this  is  analogous  to  superposing  the  two  velocity 
potentials  which,  taken  individually,  would  correctly  represent  a 
single  bubble. 

^  H.  G.  Snay,  "Underwater  Explosion  Phenomena;  The  Parameters  of 
Migrating  Bubbles,"  NAVORD  4125,  (1962) 

Compare  with  page  36  of  reference  3 
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To  first  order  we  fix  the  location  of  the  two  sources  at  the 
Initial  separation  distance  of  the  two  charges.  The  justification 
of  this  approach  Is  provided  by  the  previous  study^  of  the  pulsation 
and  migration  of  two  Interacting  spherical  bubbles.  It  was  found 
that  the  centers  of  the  Interacting  bubbles  remain  almost  stationary* 
at  first,  but  strongly  move  toward  each  other  later  as  the  bubble 
contracts.  This  assumption  Is  further  substantiated  by  the  good 
agreement  between  calculations  based  on  Equation  (2.^1)  and  calcula¬ 
tions  with  the  MACYL  Code  (See  section  3,2). 

2.2  Equations  for  the  Interface.  Since  a  bubble-water  Interface 
consists  of  the  same  particles  throughout  Its  motion**,  the  position 
of  a  particular  point  on  the  bubble  Interface  can  be  determined  by 
Integration  of  the  particle  velocities  over  time. 

The  particle  velocities  for  any  point  In  the  flow  field  are 
given  by  the  gradient  of  the  potential  given  by  Equation  (2.il), 


(2.5a) 


(2.5b) 


In  fact,  the  two  equal  spherical  bubbles  slightly  repel  one 
another,  but  the  total  distances  moved  are  small. 

Vaporization  or  boiling  of  the  ambient  water  Into  the  bubble 
due  to  passage  of  the  shock  wave  Is  ignored. 
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where  u  and  v  are  the  particle  velocities  in  the  x  and  y  directions, 
respectively,  as  shown  in  Figure  1. 

Thus,  at  any  time  t,  the  coordinates  x  and  y  of  a  point  on  the 
bubble  interface  are  given  by 


and 


x(t)  »  x(to)  +  ^  u  dt 


t. 


(2 . 6a) 


=  x(to) 


t  f  <3(t)  /jL  +  (as±Bi\ 

I  <  ) 


dt 


y(t)  =  y(to)  +  y*  v  dt 

to 


=  y(to) 


+  f  Q(t)y/ dt, 

I  ^7 


(2.6b) 


where  x(to)  and  y(to)  are  the  coordinates  of  the  interface  particle 
at  time  to.  It  is  seen,  then,  that  specifying  Q(t)  is  sufficient 
to  find  the  motion  of  an  Interface  particle.  By  tracing  the  trajec¬ 
tories  of  enough  interface  particles,  it  is  possible  to  numerically 
determine  the  motion  of  the  whole  interface  and  obtain  the  resulting 
flow  field. 

2.3  Determination  of  Q(t).  The  function  Q(t)  may  be  determined 
from  the  equation  for  the  conservation  of  energy. 

The  kinetic  energy  of  the  water  is  given  by  integration  of  the 
square  of  the  velocities  over  the  volume  of  the  field. 
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flow  field 
volume 


■  K/j  (u)'-  (»)■}..  ■ 

flow  field 
volume 


(2.7) 


By  use  of  Green’s  formula,  this  volume  integral  may  be  reduced  to 
a  surface  integral  enclosing  the  flow  field  volume.  As  there  is  no 
flow  at  infinite  distance,  the  surface  Integral  vanishes  at  infinity. 
This  means  that  the  surface  Integral  needs  to  be  taken  over  the 
bubble  Interface  o  only.  Thus,  the  kinetic  energy  in  the  flow 
field  becomes 


T  = 


(2.8) 


bubble 

Interface 

where  3(j)/3n,  the  partial  derivative  of  in  the  direction  of  the 
outward  normal  to  the  Interface,  is  the  normal  velocity  and  Po  is 
the  density  of  the  ambient  water.  Equation  (2,8)  can  be  expressed 
in  terms  of  Q(t)  by  use  of  Equations  (2.5): 


loItT  ^  5TtT  S' )  •  " 

bubble 

Interface 

'  / 

*■  . 

(2.9a) 


/S  A 

Where  x  and  y  are  the  unit  vectors  in  the  x  and  y  directions  and  n 
is  the  outward  pointing  unit  vector  normal  to  the  bubble  Interface. 


f 

t 

t 
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The  expression  In  the  square  brackets  Is  the  dot  product.  For  a 
single  spherical  bubble  the  kinetic  energy  becomes 

T  =  2TrpoA^A=^  =  27rpoQVA.  (2.9b) 

The  potential  energy,  ,  of  one  of  the  two  bubbles  Is  given  by 

^  =  PqV  ,  (2.10) 

where  Pq  is  the  hydrostatic  pressure  In  the  water  and  V  Is  the 
volume  of  the  bubble. 

The  gaseous  contents  of  the  bubbles  are  taken  as  Ideal  gases 
with  a  constant  adiabatic  exponent  y ,  an  assumption  carried  over 
from  the  classic  bubble  theory.  For  each  bubble,  the  Internal  energy 
of  the  gaseous  bubble  contents,  G,  Is  given  by 

“  = '“r.  (v)  "  .  *  C2.ll) 

where  and  are  arbitrary  reference  magnitudes,  for  Instance 
the  Initial  pressure  and  the  corresponding  Initial  volume. 

The  total  energy  of  the  mechanical  system  comprising  the  two 
bubbles  and  the  surrounding  water,  2E,  Is  conserved.  Thus,  the 
sum  of  the  kinetic,  potential,  and  Internal  energy  Is  constant 
throughout  the  first  pulsation.  The  magnitude  E  Is  called  the 
bubble  energy.  For  a  single  bubble.  It  Is  given  by 

T  +  +  G  =  E  (2.12a) 

or 

*  This  equation  holds  only  for  two  bubbles  produced  by  the 
same  explosive  and  charge  weight.  If  the  type  of  explosive  or 
the  charge  weight  is  not  the  same,  a  different  expression  for  G 
Is  required  for  each  bubble.  When  such  bubbles  coalesce,  their 
pressures  may  not  be  the  same  and  a  pressure  equalization  takes 
place  at  the  moment  of  their  combination. 
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T  +  PoV  + 


E. 


(2.12b) 


The  desired  expression  for  the  source  strength  Q(t)  of  one  bubble 
may  be  obtained  from  Equations  (2.12)  and  (2.9).  At  the  Instant  t  we 
have 


Q(t)  = 


^(E-n-o)  /  i(,„y J  " 


bubble 

Interface 


do. 


(2.13) 


where  the  positive  square  root  is  used  as  the  bubbles  are  expanding. 
Equations  (2.12)  and  (2.13)  may  be  somewhat  simplified  by  using 
non-dimensional  coordinates  due  to  Friedman^.  The  transform  to 
dimensionless  variables  is  based  upon  the  bubble  energy  E.  The 
length  scale  L  and  the  time  scale  C  are  given  as 


L  =  (3E/il7TPo)^^'^ 

=  1.7289  (eW/Z)^/3 


(2.1i|a) 


C  =  L(3Po/2Po)^^^ 

(2.l4b) 

=  0.37331  (eW)^^^Z"^^^  , 

where  P©  Is  the  hydrostatic  pressure,  W  Is  the  weight  of  one  of  the 
explosive  charges,  Z,  the  hydrostatic  pressure  expressed  in  feet  of 
ambient  water,  and  e,  the  specific  bubble  energy  in  cal/gra.  Intro¬ 
ducing  the  dimensionless  variables 
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t  =  t/C 
A  =  A/L 
r  -  r/L 


a  = 
V  = 
Q  = 
T  = 


V/  ijL* 
Q/L*(L/C) 


T/2irpoL* 


we  obtain  from  Equation  (2.12) 


T  +  V  +  kV 


(1-Y). 


1, 


(2.1ilc) 


(2.15) 


where  k  Is  the 


dimensionless  bubble 
/ 


parameter  defined  as 

Ty^I- 


(2.16) 


Now  Equation  (2.13)  for  dimensionless  source  strength  Q  becomes 


Q(t) 


/ 


[d-v-k  V 


(l-T) 


bubble 
Interface 


do 


(2.17) 


Selection  of  Bubble  Parauneters  and  Initial  Conditions.  In 

order  to  find  Q(t)  and  the  position  of  the  points  of  the  Interface 

by  means  of  Equation  (2.6),  It  Is  necessary  to  specify  y  for  the 

bubble  contents,  the  Initial  Interface  shape,  as  well  as  P  and  V  . 

r  r 

These  parameters  have  been  the  subject  of  considerable  study.^*^ 
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The  values  of  the  adiabatic  exponent  y  for  the  gaseous  contents  of 
the  bubble  are  known  for  several  types  of  high  explosives.  If  the 
dimensionless  notation  is  usedj  the  determination  of  and  is 
tantamount  to  determining  k  (provided,  of  course,  that  y  and  the 
depth  are  given).  On  the  basis  of  experimental  pressure  measure- 
ments  made  using  single  TNT  explosive  charges,  Arons'  gives  the 
following  expression  for  k, 


k  =  0.0552 


(2.18) 


Strictly  speaking,  this  relation  should  not  be  used  if  the  depth  of 
explosion  is  greater  than  1000  feet  but  this  is  not  a  serious  limi¬ 
tation  for  our  consideration. 

Initially,  the  two  bubbles  are  assumed  to  be  spherical.  Their 
radius  can  be  determined  from  Equation  (2.15) •  Just  after  detonation 
the  dimensionless  kinetic  energy  in  the  water  T  is  very  small. 
Further,  the  dimensionless  volume  V  is  also  small.  Equation  (2.15) 
becomes 

k  =  1 

(2.19) 

or 

Vo  = 


The  dimensionless  initial  radius,  Aq ,  may  thus  be  obtained  from  Voas 


Ao  =  (kl/^Y-l))l/3  (2.20) 


^  A.  B.  Arons,  "Secondary  Pressure  Pulses  Due  to  Gas  Globe  Oscilla¬ 
tions  in  Underwater  Explosions  II.  Selection  of  Adiabatic  Parameters 
in  the  Theory  of  Oscillation,"  Journal  of  the  Acoustical  Society  of 
America,  2^,  (19^8).  Reprinted  in  Underwater  Explosion  Research, 
Volume  II,  The  Gas  Globe,  Office  of^Naval  Research,  (1950 ) 
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which  Is  analogous  to  a  similar  expression  In  the  classic  incompress¬ 
ible  theory  for  a  single  bubble*. 

The  time,  to,  at  which  the  bubbles  are  spherical  Is  taken  as 
to=0,  l.e.,  the  Instant  of  detonation.  Strictly  speaking,  any  Incom¬ 
pressible  representation  of  explosion  phenomena  Is  not  valid  for 
times  very  shortly  after  detonation.  For  these  short  times, 
phenomena  which  are  due  to  the  compressibility  of  the  surrounding 
water  (l.e.,  the  shock  wave  formation)  are  dominant.  An  Incompress¬ 
ible  theory  cannot  possibly  be  in  accord  with  reality  at  these  times. 
It  Is  not  surprising  then  that  the  radius  given  by  Equation  (2,20) 
does  not  correspond  to  the  radius  of  the  explosive  charge;  Equation 
(2.20)  gives  a  radius  about  70^  too  large.  When  the  values  for 
Ao  and  k  from  Equations  (2.18)  and  (2.20)  are  used  to  represent  a 
single  explosion  bubble,  the  results  for  the  maximum  radius  and 
period  agree  with  experimental  measurements  to  within  a  few  percent. 
These  calculations  also  yield  histories  of  the  bubble  radius  and  the 
pressure  in  the  water  which  agree  favorably  with  experimental  re¬ 
cords.  There  Is  every  reason  to  believe  these  initial  conditions 
are  also  applicable  to  explosion  bubbles  from  two  separated  charges. 

The  Computer  Program.  The  developments  discussed  In  Sections 
2.2,  2.3,  and  2.4  were  written  into  a  computer  code  which  calculates 
the  expansion  of  two  non-spherlcal  bubbles.  The  time  integrations 
for  the  coordinates  of  the  bubble  interface,  Equation  (2.6),  were 
done  numerically.  Several  interface  particle  trajectories  must  be 
calculated  so  that  at  each  time  step  the  entire  interface  shape  can 

be  determined  by  Interpolation.  The  bubbles  are  axlsymmetrlc  about 

K - - - 

See  Equation  (13)  of  reference  7. 
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the  x-axls  of  Figure  1.  This  means  that  the  particle  trajectories 
in  only  one  plane  containing  the  axis  of  symmetry  (e.g.,  the  x  -  y 
plane  of  Figure  ,1)  need  be  considered.  The  interpolations  were  done 
using  sections  generated  by  rotation  of  segments  of  parabolas  in  the 
X  -  y  plane  about  the  x-axls  near  the  y-axls.  The  directrix  of  the 
parabolas  selected  was  along  the  y-axls.  Near  the  x-axis  and  the 
median  plane  (l.e.,  the  plane  equidistant  from  both  charges)  the 
directrix  was  taken  along  the  x-axls.  Thus,  at  each  time  step  the 
volume  of  the  bubbles  and  the  surface  integration  of  Equation  (2.17) 
is  then  used  to  find  Q(t)  at  this  time.  This  value  of  Q  is  then 
used  in  the  quadratures  of  Equation  (2.6)  to  find  the  new  position 
of  interface  particles,  etc.  The  time  integrations  were  made  using 
the  PNOL3  computer  program  described  in  a  report  by  Ferguson  and 
Orlow® . 

Description  of  the  MACYL  Calculations.  The  MACYL  computer  pro- 

9 

gram  is  a  technique  for  solving  incompressible,  viscous,  two- 
dimensional,  axially  symmetric,  initial-value  fluid  flow  problems 
involving  zero,  one,  or  two  free  surfaces.  The  incompressible 
Navler-Stokes  equations  and  the  continuity  equation.  Equation  (2.21), 
are  solved  by  finite  difference  methods. 


11^  +  (V.grad)V  =  -  grad  p  +  J-  V 

M  P 


(2.21) 


R.  E.  Ferguson  and  T.  A.  Orlow,  ”FN0L3,  A  Computer  Program  to 
Solve  Ordinary  Differential  Equations,”  NOLTK  71-2,  (1971) 

M  J*  P*  Shannon,  B.  J.  Daly,  The  MAC 

Method— A  Computing  Technique  for  Solving  Viscous,  Incompressible, 
Transient  Fluid-Plow  Problems  Involving  Free  Surfaces,"  Los  Alamos 
Scientific  Laboratory,  LA-3425,  (1966) 
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+  p  dlv  V  +  Vgrad  p  =  0 

where  n  Is  the  dynamic  viscosity.  For  details  of  the  computing 
method  see  references  9  and  10.  A  sample  of  the  MACYL  output  Is 
shown  In  Figure  2 .  This  figure  shows  the  shape  at  the  maximum 
expansion  of  a  single  bubble.  The  slight  deviation  from  the 
spherical  shape  seen  is  due  to  the  effect  of  gravity.  The  "scatter" 
of  the  points  Is  typical  of  such  calculations.  It  has  been  Improved 
considerably  since  these  calculations  were  made. 

3 .  NUMERICAL  CALCULATIONS 

3.1  Present  Method.  Calculations  were  made  assuming  two  347-lb 
TNT  charges  with  separations  of  10,  30,  and  100  ft  exploding  In  a 
depth  of  70  ft.  The  shapes  of  the  Interacting  bubbles  are  Illus¬ 
trated  In  Figures  3,  4,  and  5.  Typically,  20  interface  particle 
paths  were  considered.  The  figures  show  the  trajectories  of  these 
particles.  This  number  of  Interface  paths  gives  agreement  within 
1%  for  the  maximum  radius  and  the  half  period  of  spherical  bubbles. 

To  obtain  this  check,  the  bubble  pulsation  was  calculated  by  means 
of  the  classical  theory  (Incompressible  fluid)  for  charge  weights 
of  347-lb  TNT  and  694-lb  TNT. 

The  Initial  conditions  and  the  Important  magnitudes  associated 
with  these  calculations  were  as  follows: 

J.  W.  Pritchett,  "MACYL — A  Two-Dimensional  Cylindrical  Coordinate 
Incompressible  Hydrodynamic  Code,"  Naval  Radiological  Defense  Labora¬ 
tory  USNRDL-LR-67-97,  (1967) 
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Depth  of  Explosion 
Isentroplc  Exponent 
Bubble  Parameter 

Bubble  Energy 
Length  Scale  Factor 
W  =  3^  7-lb  TNT 
W  =  69^-lb  TNT 
Time  Scale  Factor 
W  =  347-lb  TNT 
N  =  694-lb  TNT 
Initial  Radius 

Dimensionless 
W  =  347-lb  TNT 
W  =  694-lb  TNT 
Initial  Pressure 


D  =  70  ft 
Y  =  1.2415 
k  =  0.0552 
=  0.16906 

e  =  490  cal/gr 

L  =  1.7289  (eW/Z)^/3 
=  20.433  ft 
=  25.744  ft 
0  =  0.37331 

=  0.4347  sec 

=  0.5477  sec 


Ao  -  Aq/L 
Aq 
Ao 
P(Ao) 


=  =  0.08599 

=  1.7570  ft 
=  2.2137  ft 

=  17,426.6  psl 


Single  Spherical  Charge: 

Maximum  Radius  Half  Period* 

W  =  347-lb  TNT  19.096  ft  0.327  sec 

W  =  694-lb  TNT  24.06  ft  0.412  sec 

3.2  MACYL  Calculations.  The  results  of  calculations  using  the 
method  in  this  report  were  compared  with  the  calculations  using  the 
MACYL  Code.  Figure  6  Illustrates  the  agreement  between  these  two 
techniques.  This  figure  is  drawn  for  a  charge  separation  of  10  ft 

and  for  two  347-lb  TNT  charges.  The  method  starts  with  the  same 

TF— - — — - 

Half  period  is  the  time  of  the  bubble  maximum. 
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Inll/ldrl  conditions  as  those  used  in  our  calculations.  Agreement 
between  the  two  methods  for  the  expansion  phase  of  the  double 
bubble  is  seen  to  be  within  10%. 

4 .  DAMAGE  EFFECTIVENESS 

4-1  Bubble  Parameters  Which  Cause  the  Whipping  of  Ships.  It  was 
originally  believed  that  a  double  charge  arrangement  would  produce 
greater  damage  to  naval  targets  than  a  single  charge  of  equal  total 
weight.  The  analysis  presented  here  has  shown  that  this  is  not  the 
case.  In  this  paragraph,  we  present  a  summary  of  the  considerations 
which  originally  led  to  this  assumption.  Subsequently,  in 
paragraph  4.3,  It  will  be  shown  that  the  results  of  the  present 
study  on  the  behavior  of  double  bubbles  invalidates  the  basis  for 
this  assumption. 

Since  we  are  dealing  with  bubble  phenomena,  an  underbottom 
attack  against  surface  ships  is  considered,  as  this  Is  the  most 
Important  type  of  damage  done  by  the  underwater  explosion  bubble. 

In  this  case,  the  explosion  is  set  off  at  some  distance  below  the 
bottom  of  the  ship.  The  expanding  bubble  excites  violent  trans¬ 
verse  flexural  vibrations  of  the  entire  body  of  the  ship,  known 
as  whipping.  These  vibrations  almost  always  occur  In  the  vertical 
plane,  l.e.,  they  are  up-and-down  motions. 

The  ship  may  be  visualized  as  a  floating  beam  (a  girder)  of 
varying  cross  section.  The  dangerous  effects  of  whipping  arise  from 
the  different  acceleration  of  the  sections  of  the  ship.  The  section 
above  the  explosion  Is  accelerated  upward  most  strongly.  The 
acceleration  decreases  with  Increasing  longitudinal  distance  from 
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the  attacked  section.  This  uneven  acceleration  produces  a  bending 
of  the  ship  and  Induces  strong  flexural  oscillations^^  which  can 
severely  Impair  the  girder  strength  of  the  ship  and  may  even  sever 
the  ship  Into  two  parts. 

An  explosion  beneath  the  ship  produces  rapid  motions  of  the 
surrounding  water.  This  flow  exerts  hydrodynamic  forces  on  the 
ship  which  are  proportional  to  the  average  vertical  acceleration 
il  of  the  water  displaced  by  the  ship.  These  forces  produce  an 
upward  acceleration  y  of  the  ship  as  well  as  elastic  or  plastic 
restraining  forces  In  Its  structure. 

Consider  sections  of  the  ship,  each  bounded  by  planes  perpendi¬ 
cular  to  the  ship’s  axis.  For  the  1— th  section  we  have: 


(il.l) 


where 


\x^  =  average  upward  acceleration  of  the  water  which,  in  the 
absence  of  the  ship,  would  occupy  the  space  displaced 
by  the  1-th  section 

y^  =  bodily  upward  acceleration  of  the  1-th  section 
F^  =  total  of  the  restraining  forces  acting  on  both  faces  of 
the  1-th  section 

mpi  =  mass  of  water  displaced  by  the  1-th  section 
"'wi  “  ^<^ded  mass  of  the  1-th  section  =  mass  of  water  which 
appears  to  move  with  the  section.  depends  on  the 


Elastic  Theory  of  the  Explosion  Induced  Whipping 
Motion  of  Ships  and  Submarines"  NCRE/R.532,  (1968) 
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underwater  shape  of  the  section 
=  actual  mass  of  the  1-th  section  =  mass  of  hull,  decks, 

bulkheads,  stiffeners,  machinery,  etc.  situated  In  the 
section . 

Note  that  commonly 

®Di  ^  >  i^.2) 

whereas  according  to  Archimedes'  principle 

N  N 

S  ^Di  “  Xrf  ^1  ’  (4*3) 

i=l  1=1 

where  N  is  the  total  number  of  sections. 

The  average  upward  acceleration  u  due  to  the  expansion  of  a 
sphei’lcal  bubble  of  radius  A(t)  Is  given  by 

1 

where 

djL  =  depth  of  bubble  center  below  the  center  of  mj^^ 

"*  slant  distance  between  these  centers  . 

Equation  (4.4)  holds  for  the  case  where  the  bubble  center  Is  on 
vertical  plane  of  symmetry.  The  general  case  can  be  readily  written 
down  but  Is  not  of  interest  here. 

The  magnitude  A^A  is  the  source  strength  Q(t)  discussed  in 
section  2  of  this  paper.  It  may  thus  appear  that  QCt),  or  specifi- 

t 

cally  Q(t),  Is  the  physical  magnitude  responsible  for  the  excitation 
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of  whipping.  Here  further  details  of  the  whipping  process  must  be 
considered  before  drawing  a  conclusion. 

Equation  (^.^)  for  the  upward  acceleration  u  Is  based  pn  the 
relations  of  Incompressible  fluid  motion.  As  stated  In  paragraph  2.3 
the  assumption  of  an  Incompressible  medium  is  adequate  for  the 
description  bubble  phenomena  except  for  the  early  times  following 
the  instant  of  detonation.  This  exception  holds  for  the  whipping 
process  to  an  even  greater  degree.  At  these  early  times,  the  target 
Is  subjected  to  the  Impact  of  the  shockwave  and  Its  afterflow.  At 
this  stage  Equations  (4.1)  and  (4.4)  lose  their  physical  meaning 
and  hence  their  validity. 

The  shockwave  produces  elastic  or  plastic  deflections  of  the 
hull  of  the  ship.  The  energy  thus  transmitted  to  the  plating  Is 
either  expended  in  the  process  of  hull  rupture  or,  If  no  rupture 
occurs.  Is  transmitted  to  the  side  of  the  hull  and  the  transverse 
bulkheads.  Equalization  of  the  upward  momentum  Imparted  by  the 
shockwave  into  the  members  of  the  ship’s  structure  ultimately  yields 
upward  motions  of  the  sections  of  the  ship  that  are  analogous  to 
whipping  motions.  These  motions  are  enhanced  by  the  afterflow  of 
the  shockwave,  i.e.,  the  flow  caused  by  the  expanding  bubble.  An 
analytic  description  of  these  early  processes  Is  difficult  and  has 
not  yet  been  accomplished. 

Whipping  theory  utilizes  an  approximation  for  these  early  phases 
which,  judging  from  the  good  agreement  between  experimental  and 
calculated  results,  seems  to  describe  the  actual  phenomena  very  well. 
It  Is  assumed  that  initially  a  velocity  Is  Imparted  to  the  ship  that 
corresponds  to  the  maximum  rate  of  the  volume  expansion  of  the  bubble, 
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hence 


yi(t=o)  =  I  /aU\ 

1  ”1  +  rowi  V  / 


max 


whereby 


y^(t=0)  =  0. 


(4.6) 


Here  yj^(t)  is  the  deflection  of  the  axis  of  the  ship  at  the 
center  of  the  i-th  section.  The  magnitude  is  the  maximum 

\  J  luaX 

source  strength  Q 

max 

It  appears  from  Equation  (4.5)  that  every  section  of  the  ship 

is  decelerated  for  t  >  0.  This  is  because  yx(t=0)  and  Q  are 

i  max 

maximum  values  and,  hence  y^(t)  and  Q(t)  must  subsequently  decrease. 

For  the  same  reason  y^^  in  Equation  (4.1)  as  well  as  Q  are  negative 

for  t  >  0.  It  is  thus  seen  that  the  basic  equation  of  the  whipping 

theory  (4.1)  considers  the  retardation  of  the  whipping  velocity  from 

its  initial  "kick-off  velocity",  ^^^(t^O).  Although  this  retardation 

is  Important  for  details  of  the  whipping  process,  it  appears  that 

Qjjjax  parameter  of  predominant  importance  for  whipping  damage. 

To  summarize,  the  maximum  source  strength  Q  is  the  primary 

max 

measure  of  whipping  damage  which  an  expanding  explosion  bubble  can 
inflict  to  a  ship.  The  negative  values  of  Q  that  occur  after  Q 

max 

are  parameters  of  secondary  significance.  The  smaller  |q|,  the  less 
decelerated  is  the  whipping  motion  and  the  larger  is  the  damage 
under  otherwise  equivalent  conditions. 

4.2  Damage  Potential  of  Double  Bubbles.  Figure  7  Illustrates  the 
superposition  of  the  velocities  produced  by  the  expansion  of  two 
separate  bubbles.  As  stated  above,  the  velocity  components  from 
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each  bubble  are  proportional  to  the  source  strength  of  each  bubble 
j  =  1,2.  To  obtain  a  simple  representation  of  Q  we  set 


and 


whence 


and 


A(t)  = 

(^.6) 

t  = 

(^.7) 

Act)  =  1  z^/2p^(T), 

(^.8) 

A^A  = 

(^.9) 

(^.10) 

Here 

J  =  Radius  Coefficient 
K  =  Period  Coefficient 
W  =  Charge  Weight 
Z  =  Hydrostatic  Head. 

The  functions  ....  p^  are  virtually  independent  of  charge  weight, 
depth,  or  the  type  of  the  explosive,  but  depend  on  the  dimensionless 
time  T  only.  The  trend  of  these  functions  is  not  Important  for  our 
study. 

According  to  Figure  7,  the  resultant  upward  velocity  u  from 
two  separated  bubbles  generated  by  the  simultaneous  explosion  of 
charges  having  the  weights  and  is 
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_  sln0  „-l/6 


(^.11) 


or  for  equal  charge  weights  with  W  =  =  W2  and  with  R  =  D/slnG 


u 


K 


3 


X  2. 


(il.lla) 


The  magnitudes  0,  D,  and  R  are  shown  in  Figure  7*  If  the  two  charges 
are  combined  into  one  single  charge,  we  obtain 


u  =  P3(t)  (4.12) 

or  If  and  are  equal 

u  =  f  2^/^  /  .  (4.12a) 

Hence  the  resultant  velocity  u  from  two  equal  separated  bubbles  is 
sln®0x2/2  times  larger  than  that  from  a  lumped  charge  of  equal 
total  weight  measured  at  the  same  distance.  This  corresponds  to  a 
26%  Increase  If  0  Is  close  to  90°. 

For  the  acceleration  u,  one  finds 


■ 

1 

u  =  — 

sln0  22/3 
R^ 

+  F,  (t 

• 

_ 

or  for  equal  charge  weights 
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u 


K» 


X 


2. 


For  a  lumped  charge  we  have 


(^.13a) 


u 


K* 


r2/3 


Fi.(t) 


(^.14) 


or  for  equal  charges  and  W2 

u  -  ^  ?/,(!)  X  2^^^  .  (4.14a) 

D* 

Here  the  acceleration  u  from  two  equal  separated  bubbles  is 
sin^e  X  2/2^^^  =  sln^O  x  1.59  times  larger  than  that  from  a  lumped 
charge  of  total  equal  weight. 

This  possible  enhancement  of  almost  60J5  originally  stimu¬ 
lated  the  study  described  in  this  paper.  It  was  quickly  recognized 

that  not  the  acceleration  u  but  the  maximum  velocity  u  must  be 

max 

used  as  a  criterion  for  the  effectiveness.  This  decreased  the  maxi¬ 
mum  enhancement  to  26X,  which  would  still  be  a  welcome  gain. 

This  enhancement  is  obtained  if  the  two  bubbles  move  independ¬ 
ently  and  do  not  interact.  If  interaction  takes  place  a  reduction 
of  the  enhancement  must  be  expected  since  the  greatest  interaction 
corresponds  to  a  combination  into  a  single  bubble.  No  enhancement 
occurs  in  this  case. 

^•3  Effectiveness  of  Interacting  Bubbles.  The  interaction  of  two 
bubbles  is  the  subject  of  our  paper.  Figure  8  shows  the  source 
strength  Q  for  one  of  the  two  bubbles  as  a  function  of  time.  The 

charge  arrangements  consist  of  two  347-lb  spherical  TNT  charges 
at  10- ,  30- ,  and  100-ft  separation.  Also  shown  is  the  source 
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strength  for  a  non-interacting  bubble  (infinite  separation).  Since 
the  charge  radius  of  a  3^7-lb  charge  Is  roughly  one  foot,  the  curves 
may  also  be  used  for  other  charge  weights  by  expressing  the  separa¬ 
tion  In  charge  radii.  More  accurately,  the  curves  hold  for  10,6, 

31. 7 j  and  105.7  TNT  charge  radii. 

Figure  8  Illustrates  that  the  Interaction  of  the  two  bubbles 
reduces  the  maximum  source  strength,  but  Increases  the  half  period. 
Even  at  the  large  separation  of  100  ft,  a  substantial  Interaction  Is 
noticeable.  For  a  10-ft  separation,  the  Interaction  has  Increased  so 
much  that  the  maximum  source  strength  is  reduced  to  that  of  a  lumped 
charge,  shown  by  the  dashed  curve  In  Figure  8. 

The  dashed  curve  In  Figure  8  Is  obtained  as  follows:  The  solid 
curves  show  the  source  strength  Q  for  one  of  the  two  interacting 
bubbles.  Being  symmetric  the  other  bubble  produces  the  same  Q,  so 
that  the  total  strength  is  2Q.  If  the  separation  B  is  reduced  to 
zero  a  lumped  charge  of  694-lb  TNT  is  obtained.  Its  total  source 
strength  Is  2Q  in  our  notation.  To  be  consistent  with  the  curves  of 
Figure  8,  we  have  shown  ^(694  lb)  for  the  source  strength  of  the 
lumped  charge. 

The  source  strength  Q  is  proportional  to  the  square  of  the 
length  scale,  which  in  turn  is  proportional  to  the  cube  root  of  the 
bubble  energy  or  the  charge  weight.  This  is  seen  from  the  definition 
of  the  dimensionless  source  strength  Q  given  in  Equation  (2.l4c). 
(Note  that  the  time  scale  factor  C  is  proportional  to  the  length 
scale  factor  L.  Hence,  C/L  which  occurs  in  the  definition  of  Q  is 
independent  of  the  charge  weight.)  It  follows  that 
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^(69^1  lb)  =  I  1.5874Q(3^7  lb)  -  0.7937Q(3^7  lb). 

Since  the  dashed  curve  of  Figure  8  refers  to  a  single  explosion,  it 
can  be  compared  with  the  curve  for  infinite  separation  which  corre¬ 
sponds  to  Q(3^7  lb).  Therefore,  the  ordinates  of  the  dashed  curve 
are  79*37?  of  those  of  Q(3^7  lb),  whereas  its  time  scale  is  2^^^  ■ 
1.25992  larger. 

If  one  adopts  the  maximum  source  strength  as  a  measure  of  the 
severity  of  whipping  damage,  it  is  seen  that  no  advantage  in  damage 
effectiveness  is  apparent  for  a  double  charge  of  ten-feet  separation: 
For  sln0  ■  1,  the  resultant  maximum  upward  velocity  is  Just  equal  to 
that  produced  by  a  lumped  charge.  One  may  even  argue  that  separated 
charges  are  inferior  because  in  practical  applications  sinG  is  less 
than  unity.  This  argument  holds  particularly  true  for  the  30-ft  or 
100-ft  separation.  If  the  normal  standoff  from  the  ship's  bottom  is 
30  ft,  we  have  0  =  63.^°  and  sln*0  =  0.7155  for  30-ft  charge  separa¬ 
tion  and  0  =  31.0°  and  sin*0  »  O.1362  for  100-ft  separation.  In  both 

cases  the  term  sln’G  more  than  compensates  for  the  increase  of  Q 

^ax 

shown  in  Figure  8,  which  is  8?  for  30-ft  and  19?  for  100-ft  separation. 

In  the  above  considerations  the  resultant  velocity  on  the 
perpendicular  bisector  is  considered.  In  actuality,  the  velocity 
distribution  along  the  ship  must  be  taken  into  account.  Here  multiple 
charges  have  the  advantage  that  the  velocity  decreases  less  from  its 
maximum  value.  Hence,  the  velocity  profile  along  the  bottom  of  the 
ship  is  fuller  for  multiple  charges  than  for  a  single  charge.  But  a 
single  charge  usually  produces  a  higher  maximum  velocity  on  the 
normal  bisector. 
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This  effect  can  only  be  accounted  for  by  means  of  whipping 
calculations.  Such  a  study  has  been  completed  at  the  Naval  Ordnance 
Laboratory  and  It  has  confirmed  the  conclusion  of  the  present  study 
that  no  advantages  can  be  expected  from  multiple  charges  in  doing 

whipping  damage  to  ships.  These  calculations  will  be  described  in  a 
forthcoming  paper, 

5 .  CONCLUSIONS 

V/e  have  demonstrated  that  potential  flow  theory  using  a  single 
stationary  point  source  of  variable  strength  Is  sufficient  for  the 
calculation  of  the  shape  and  the  properties  of  two  Interacting 
bubbles  during  the  expansion  phase.  The  results  of  our  simple 
calculations  agree  within  10^  with  those  obtained  by  finite  differ¬ 
ence  methods  using  the  MACYL  computer  technique.  The  justification 
for  assuming  a  stationary  (not  translating)  point  source  is  given  by 
the  results  of  our  previous  study^,  which  show  that  the  centers  of 
spherical  bubbles  expanding  in  unison  remain  essentially  stationary. 
Upon  contraction  (not  considered  here)  the  centers  move  rapidly 
toward  each  other. 

The  maximum  strength  of  the  source  representing  the  expanding 
bubble  Is  a  measure  for  the  degree  of  whipping  damage  done  by  the 
bubble. 

The  interaction  between  two  expanding  bubbles  reduces  the 
maximum  source  strength  even  for  large  separations  of  bubble  centers. 
As  a  consequence,  double  explosions  offer  no  advantages  in  bubble 
damage  effectiveness  against  ships. 
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R  ANDZ  IN  FEET 

FIG.  2  SAMPLE  OUTPUT  OF  A  CALCULATION  MADE  USING  THE  MACYL  COMPUTER  TECHNIQUE 
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The  Calculations  Were  Made  For  An  Explosion  Of  Two  347-lb  TNT  Charges 
10  Feet  Apart  At  A  Depth  Of  70  Feet.  Coalescence  Occurs  Between  The 
Reduced  Time  Of  0.  05  And  0. 16. 
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FIG.  4  BUBBLE  INTERFACES  AND  PARTICLE  PATHS  FOR  DOUBLE  EXPLOSION  DUBaLES 
SEPARATED  BY  31 .7  CHARGE  RADII 

Tlie  Colculatfoni  Woro  Made  For  An  Explosion  Of  Two  347-lb  Charges  30  Ft. 

Aporf  At  A  Depfh  Of  70  Ft. 


T  ■  TIME/TIME  OF  MAXIMUM  EXPANSION 


NOLTR  73-102 


li 


I 

I 


UJ 


UI  U 


QC 

< 


eo 

A 

z 

g 

o 

is; 

X 


to 

3 

o 

Q 

ec 

O 

u. 

«/> 

X 

H* 

UJ  S 

da 


Z  =  ^ 
.9  3 

ill 

Sjif 

$  ^  u 

tic 

o 

§  §!  S 

ill 


fcuj 

<<5 

Ox 

ZU 

<  *-  i 

5  ^ 

'a-  •' 

<^^ 

i  .  i 

in 
UJ  IQ 

S  i  8- 

Zji; 

2^3 
.2  5  x  » 
-■=1  » 

m  < 

«S 

5**  «- 

3  Ul 
eO  U1 

itl  S 

m 

1—  <  < 

o 


32 


NOLTR  73-102 


- POTENTIAL  SOLUTION 

- MACYL  CODE 


FIG.  6  COMPARISION  OF  CALCULATIONS  MADE  WITH  POTENTIAL  FLOW  THEORY 
WITH  CALCULATIONS  BY  THE  MACYL  CODE 
The  Curves  Show  One  Quadrant  Of  The  Bubble  Expansion  Of  Two  Simulta¬ 
neously  Initiated  347-lb  TNT  Charges  Separated  By  10  Feet,  The  Burst 
Depth  Is  70  Feet.  The  Parameter  Is  The  Time  In  Seconds.  The  Maximum 
Expansion  Is  Reached  At  A  Time  Of  0. 405  Seconds. 
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FIG.  7  EXPLOSION  GEOMETRY  BELOW  TARGET 
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10,000 


B-  100  FT 


B-  30  FT 


B=  lOFT 


CHARGE  WEIGHT  =  347  LB  TNT 
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DASHED  CURVE 
SHOWS  a^2  FOR  A 
694-LB  TNT  CHARGE 


LUMPED  CHARGE 


AND  B«  10  FT 


U.1  0.2  0.3  a- 

TIME  (SEC) 

8  SOURCE  STRENGTH  OF  TWO  INTERACTING  BUBBLES  AS  A  FUNCTION  OF  TIME 

The  Solid  Curves  Give  The  Source  Strength  Of  One  Of  The  Two  Interacting 
Bubbles.  The  Dashed  Curve  Shows  The  Corresponding  Magnitude  For  A 
Lumped  Charge.  The  Maximum  Source  Strength  Is  Considered  A  Measure  Of 
Damage  Effectiveness  For  Whipping.  The  Half  Period  Is  The  Time  at  Which 
The  Source  Strength  Vonlshes.  The  Graph  Shows  That  The  Maximum  Source 
Strength  Decreases  With  Decreasing  Charge  Separation.  For  A  10-Ft 
Separation  It  Assumes  The  Value  Corresponding  To  A  Lumped  Charge. 


